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(57) Abstract: A cardiac rhythm management system predicts when an arrhythmia will occur and in one embodiment invokes a 
therapy to prevent or reduce the consequences of the arrhythmia. A cardiac arrhythmia trigger/marker is detected from a patient, and 
based on the trigger/marker, the system estimates a probability of a cardiac arrhytmnia occurring during a predetermined future time 
interval The system provides a list of triggers/markers, for which detection values are recurrently obtained at various predetermined 
time intervals. Based on detection values and conditional probabilities associated with the triggers/markers, a probability estimate 
of a future arrhythmia is computed. An airhythmia prevention therapy is selected and activated based on the probability estimate of 
the future arrhythmia. 
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CARDIAC RHYTHM MANAGEMENT SYSTEM WITH ARRHYTHMIA 
PREDICTION AND PREVENTION 

Technical Field 

This document relates generally to cardiac rhythm management systems 
and particularly, but not by way of limitation, to a system providing prediction of 
a future arrhythmia and preventive therapy for avoiding or mitigating the 
5 predicted arrhythmia. 
Background 

The human heart normally maintains its own well-ordered intrinsic 
rhythm through generation of stimuli by pacemaker tissue that results in a wave 
of depolarization that spreads through specialized conducting tissue and then into 

10 and through the myocardium. The well-ordered propagation of electrical 
depolarizations through the heart causes coordinated contractions of the 
myocardium that results in the efficient pumping of blood. In a normally 
functioning heart, stimuli are generated under the influence of various 
physiological regulatory mechanisms to cause the heart to beat at a rate that 

1 5 maintains cardiac output at a level sufficient to meet the metabolic needs of the 
body. Abnormalities of excitable cardiac tissue, however, can lead to 
abnormalities of heart rhythm that are called arrhythmias. All arrhythmias stem 
from one of two causes: abnormalities of impulse generation or abnormalities of 
impulse propagation. Arrhythmias can cause the heart to beat too slowly 

20 (bradycardia, or a bradyarrhythmia) or too quickly (tachycardia, or a 

tachyarrhythmia), either of which may cause hemodynamic compromise or 
death. 

Drug therapy is often effective in preventing the development of 
arrhythmias and in restoring normal heart rhythms once an arrhythmia has 
25 occurred. However, drug therapy is not always effective for treating particular 
arrhythmias, and drug therapy usually causes side-effects that may be intolerable 
in certain patients. For such patients, an alternative mode of treatment is needed. 
One such alternative mode of treatment includes the use of a cardiac rhythm 
management system incorporated into an implantable device that delivers 
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therapy to the heart in the form of electrical stimuli. Such implantable devices 
include cardiac pacemakers that deliver timed sequences of low energy electrical 
stimuli, called pacing pulses, to the heart, via an intravascular leadwire or 
catheter (referred to as a lead) having one or more electrodes disposed in or 
5 about the heart. Heart contractions are initiated in response to such pacing 
pulses (referred to as capturing the heart). By properly timing the delivery of 
pacing pulses, the heart can be induced to contract in proper rhythm, greatly 
improving its efficiency as a pump. Pacemakers are often used to treat patients 
with bradycardia. Pacemakers are also capable of delivering paces to the heart in 

10 such a manner that the heart rate is slowed, a pacing mode referred to as anti- 
tachyarrhythmia pacing. 

Cardiac rhythm management systems also include 
cardioverter/defibrillators (ICD's) that are capable of delivering higher energy 
electrical stimuli to the heart. ICD's are often used to treat patients with 

15 tachyarrhythmias, that is, hearts that beat too quickly. Tachyarrhythmias can 
cause diminished blood circulation because the cardiac cycle of systole 
(contraction) and diastole (filling) can be shortened to such an extent that 
insufficient blood fills the ventricles during diastole. Besides the potential for 
such hemodynamic embarrassment, tachyarrhythmias can also degrade into even 

20 more serious arrhythmias such as fibrillation where electrical activity spreads 
through the myocardium in a disorganized fashion so that effective contraction 
does not occur. For example, in a particular type of tachyarrhythmia, referred to 
as ventricular fibrillation, the heart pumps little or no blood to the body so that 
death occurs within minutes. A defibrillator delivers a high energy electrical 

25 stimulus or shock to the heart to depolarize all of the myocardium and render it 
refractory in order to terminate arrhythmia, allowing the heart to reestablish a 
normal rhythm for the efficient pumping of blood. In addition to ICD's and 
pacemakers, cardiac rhythm management systems also include pacemaker/ICD's 
that combine the functions of pacemakers and ICD's, drug delivery devices, and 

30 any other implantable or external systems or devices for diagnosing, monitoring, 
or treating cardiac arrhythmias. 

Cardiac rhythm management systems incorporated into ICD's allow 
tachyarrhythmias to be automatically detected and treated in a matter of seconds. 
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Defibrillators are usually effective at treating tachyarrhythmias and preventing 
death, but such devices are not 100% effective at treating all tachyarrhythmias in 
all patients. As a result, some patients may still die even if the defibrillator 
delivers appropriate therapy. Also, some patients have frequent 
5 tachyarrhythmias, triggering frequent therapeutic shocks. This reduces the 
usable life of the implanted battery-powered device and increases the risk of 
therapy-induced complications. Furthermore, even if the device successfully 
treats the tachyarrhythmia, the patient may lose consciousness during the 
arrhythmia which can result in related serious or even fatal injuries (e.g., falling, 
1 0 drowning while bathing, car accident while driving, etc.). Thus, there is a need 
for a cardiac rhythm management system that predicts when an arrhythmia will 
occur and invokes a therapy to prevent or reduce the consequences of the 
arrhythmia. 
Summary 

1 5 The present invention relates to a system and method for predicting 

cardiac arrhythmias. In a particular embodiment, the system and method are 
implemented in an implantable cardiac device having one or more sensing 
channels for detecting conditioning events (e.g., marker/trigger events as defined 
below) and the capability of delivering some type of preventive arrhythmia 

20 therapy when conditions warrant it. 

In accordance with the invention, an arrhythmia is predicted by: 1) 
detecting a conditioning event statistically associated with the occurrence of an 
arrhythmia in a patient's heart; 2) computing a conditional arrhythmia 
probability for the conditioning event from past observations of instances in 

25 which the conditioning event occurs alone or together with an arrhythmia within 
a specified time period; 3) computing an estimated arrhythmia probability based 
upon the detected occurrence of the conditioning event; and 4) predicting the 
occurrence of an arrhythmia within a specified prediction time period if the 
estimated arrhythmia probability exceeds a specified threshold value. 

30 Conditioning events may be broadly classified into markers and triggers. 

A marker event corresponds to detected a physiological state that is statistically 
associated with occurrence of cardiac arrhythmias, but the causal relationship 
between the marker and the arrhythmia is not known. A conditioning event is 
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regarded as a trigger, on the other hand, if the event is thought to increase the 
risk of an arrhythmia occurring via a depolarization that serves as a source for 
the arrhythmia. Conditioning events may be detected on a beat-to-beat basis or 
over a longer time frame. Examples of conditioning events include a detected 
5 specific morphology of a waveform representing the electrical activity of the 
heart, a specific pattern of activation times of different areas of the heart as 
sensed by a plurality of electrodes, a specific sequence pattern of heartbeats with 
respect to time, a value of a measured physiological variable such as heart rate or 
blood pressure, or a statistic based upon a history of occurrences of conditioning 
10 events. 

In one embodiment, the conditional arrhythmia probability is calculated 
as a ratio of the number of observed instances in which the conditioning event is 
followed by an arrhythmia within a specified basic time period, to the total 
number of observed instances of the conditioning event. In that case, the 

1 5 estimated arrhythmia probability for an arrhythmia to occur within the specified 
basic time period after detection of the conditioning event is simply the 
calculated conditional arrhythmia probability. 

In another embodiment, the conditional arrhythmia probability CP is 
calculated by the expression: 

20 CP=l-e RT 

which assumes a Poisson probability distribution, where T is a measure of the 
specified prediction time period, and R is an estimate of the rate at which 
arrhythmias occur while the conditioning event is present. The rate R is a ratio 
of: 1) the number of instances in which the conditioning event is followed by an 

25 arrhythmia within a specified basic time period, to 2) the length of the basic time 
period multiplied by the total number of basic time periods in which the 
conditioning event is observed. The estimated arrhythmia probability for an 
arrhythmia to occur within the time T after detection of the conditioning event is 
again the conditional arrhythmia probability. Calculating the conditional 

30 arrhythmia probability in this manner allows the prediction time period T to 
differ from the length of the basic time period used to derive the conditional 
arrhythmia probability. 

4 
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In another embodiment, rather than basing the estimated arrhythmia 
probability upon the detection of a conditioning event, a rate at which the 
conditioning event occurs is detected over some period of time. The estimated 
arrhythmia probability is then calculated as the product of an estimated 
5 probability that a conditioning event will occur times the probability of an 
arrhythmia occurring within specified time period given the occurrence of the 
conditioning event (i.e. the conditional arrhythmia probability). Thus: 

estimated arrhythmia probability = (1- e* RT )( 1- e" CT ) 

10 where T is a measure of the specified prediction time period, R is an estimate of 
the rate at which arrhythmias occur while the conditioning event is present, and 
C is an estimate of the rate at which the conditioning event occurs. 

Another way of deriving a conditional arrhythmia probability, especially 
for trigger-types of conditioning events (although it can be used with any type of 

1 5 conditioning event), is to designate a particular detected trigger event as being 
responsible for causing a detected arrhythmia. Such culpability may be assigned 
based, e.g., upon the proximity in time between the trigger event and the onset of 
the arrhythmia, the magnitude of the detected trigger, or the frequency of 
occurrence of the trigger event within a specific time period prior to the onset of 

20 the arrhythmia. A conditional arrhythmia probability CP for that trigger event 
can then be calculated as a ratio of the number of instances in which the trigger 
event was deemed culpable for causing an arrhythmia, to the total number of 
occurrences of the trigger event. Also, as above, rather than basing the estimated 
arrhythmia probability upon the detection of the trigger event, a rate at which the 

25 trigger event occurs can be detected over some period of time. The estimated 
arrhythmia probability is then calculated as the product of an estimated 
probability that a trigger event will occur times the probability CP of an 
arrhythmia occurring within a specified time period T given the occurrence of 
the trigger event. Thus: 

30 estimated arrhythmia probability = CP x ( 1 - e" 01 ) 

In a preferred embodiment, the prediction of arrhythmias is based upon a 
plurality of the same or different detected conditioning events. A composite 
estimated arrhythmia probability is then computed as a combination of the 

5 
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estimated arrhythmia probabilities derived for each separately detected 
conditioning event. The separately detected conditioning events may be separate 
occurrences of the same or different conditioning events. As before, the 
composite arrhythmia probability is compared with a threshold value in order to 
5 predict the occurrence of an arrhythmia. In one embodiment, the composite 
arrhythmia probability is calculated by adding the individual estimated 
arrhythmia probabilities derived for each detected conditioning event, which thus 
assumes each individual arrhythmia probability to correspond to an independent 
event. In other embodiments, specific combinations of detected conditioning 
1 0 events are mapped in a non-linear fashion to estimated arrhythmia probabilities 
that can be added or otherwise combined with other estimated arrhythmia 
probabilities to give a composite value. In still other embodiments, the 
estimated arrhythmia probability is computed from a combination of conditional 
arrhythmia probabilities derived using different basic time periods but for the 
1 5 same prediction time period. 

The past observations of the occurrences of conditioning events and 
arrhythmias from which the conditional arrhythmia probabilities are derived can 
be taken from either population data or from data collected in real-time from a 
particular patient. In a preferred embodiment, the conditional arrhythmia 
20 probabilities are based initially upon past observations of the occurrences of 
events and arrhythmias taken from population data, and each probability is 
subsequently updated from a previous value to a present value with observations 
taken in real-time from a particular patient. In one embodiment, a conditional 
arrhythmia probability is updated only if the present value differs by a 
25 predetermined amount from the previous value. In another embodiment, the 
amount by which the present value differs from the previous value is tested for 
statistical significance before a conditional arrhythmia probability is updated. In 
another embodiment, the previous value of the conditional arrhythmia 
probability is incremented or decremented by a specific amount after a prediction 
time period in accordance with whether the arrhythmia occurred or not, 
respectively. 

In still another embodiment, the statistical association between the 
conditioning event and the occurrence of an arrhythmia is periodically 
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reevaluated using the most recent patient-specific data. If the statistical 
association (e.g., as a calculated from a chi-square test) is found to be below a 
specified value, the use of that conditional arrhythmia probability in deriving a 
composite estimated arrhythmia probability is discontinued. 
5 As aforesaid, one embodiment of the invention involves delivering 

preventive arrhythmia therapy if the estimated probability of an arrhythmia 
occurring within a specified time period (i.e., the composite estimated 
arrhythmia probability) exceeds a threshold value so that an arrhythmia can be 
predicted with some degree of certainty. Examples of such therapies capable of 

10 being delivered by an implantable device include the delivery of pharmacologic 
agents, pacing the heart in a particular mode, delivery of 
cardioversion/defibrillation shocks to the heart, or neural stimulation (e.g., 
stimulation of either the sympathetic or parasympathetic branches of the 
autonomic nervous system). Another type of therapy capable of being delivered 

15 by an implantable device in accordance with the invention (interpreting the term 
"therapy" in a broad sense) is issuance of a warning signal that an arrhythmia 
has been predicted, which warning signal may take the form of a audible signal, 
a radio-transmitted signal, or any other type of signal that would alert the patient 
or physician to the possibility of an impending arrhythmia. Such a warning 

20 signal would allow the patient to take precautionary measures and/or allow a 
treating physician to take other therapeutic steps if deemed appropriate. In 
accordance with the invention, the selection of a particular therapy to be 
delivered or not is based upon ascertaining a physiologic state of the patient and 
deciding whether or not a particular available modality of therapy is appropriate 

25 for delivery to the patient in that physiologic state. 

The modality selection process in one embodiment takes the form of a 
matrix mapping of a state vector representing a specific physiologic state (as 
determined by, e.g., the particular conditioning events used to make the 
arrhythmia prediction, the prediction time period for the estimated arrhythmia 

30 probability, the presence or not of specific conditioning events within a specified 
prior time period, or the magnitude and/or presence other detected and/or 
calculated variables) to a specific therapy or therapies considered most 
appropriate for preventing an arrhythmia in that instance (i.e., a point in "therapy 

7 
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space"). The matrix mapping is performed using a therapy matrix containing 
information relating to whether or not (and/or to what extent) a particular therapy 
modality is expected to be effective for a given physiologic state. The elements 
of the therapy matrix may thus constitute variables representing the 
5 appropriateness of a particular therapy modality given the presence of a specific 
element in the physiological state vector. 

In one particular embodiment, a prediction scheduler makes separate 
predictions for each time period in which an arrhythmia may or may not be 
predicted to occur (i.e., the prediction time period) and then makes a therapy 

10 decision using a therapy matrix specific for that prediction time period. In this 
manner, a therapy decision for a given therapy modality may be made with 
respect to the time period appropriate for that modality. For example, some 
therapy modalities can be expected to be effective in a short time period and are 
thus suitable for preventing an arrhythmia with a short prediction time period, 

1 5 while others would not be expected to be effective until after a longer time 

interval and therefore would only be suitable if the prediction time period were 
commensurate. In another embodiment, therapy decisions and predictions are 
made at time intervals without regard to the arrhythmia prediction periods, and 
the therapy matrix takes into account the prediction time period for the estimated 

20 arrhythmia probability. For example, the prediction time period may be 
incorporated into the physiologic state vector, and the therapy matrix then 
contains information related to the appropriateness of each available therapy 
modality for a given prediction time period. 

Rather than estimating a probability for the occurrence of any arrhythmia, 

25 separate estimated arrhythmia probabilities can be computed for different types 
of arrhythmias in accordance with the invention. By computing such separate 
probabilities, a more informed decision as to what mode of therapeutic 
intervention to employ may be made. The most beneficial set of separate 
estimated arrhythmia probabilities for any particular patient be expected to vary. 

30 Accordingly, selection of the arrhythmia types for which separate estimated 

arrhythmia probabilities are computed can be done dynamically by tabulating the 
number of detected occurrences of each specific type of arrhythmia and 
computing separate estimated arrhythmia probabilities for those arrhythmia types 
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occurring most frequently in the patient. Separate estimated arrhythmia 
probabilities can also be computed for different types of trigger events. Such 
trigger events can be expected to be patient-specific, and the particular mix of 
trigger events for which separate arrhythmia probabilities are calculated can be 
5 selected dynamically as described above. 
Rrifif Description of the Drawings 

In the drawings, like numerals describe substantially similar components 
throughout the several views. Like numerals having different letter suffixes 
represent different instances of substantially similar components. 
10 Figure 1 is a schematic drawing illustrating generally one embodiment of 

portions of a cardiac rhythm management system and an environment in which it 
is used. 

Figure 2 is a schematic drawing illustrating generally one embodiment of 
a cardiac rhythm management device coupled by leads to a heart. 
1 5 Figure 3 is a schematic diagram illustrating generally one embodiment of 

portions of a cardiac rhythm management device, which is coupled to a heart 
and/or other portions of the patient's body. 

Figure 4A is a block diagram illustrating generally one conceptual 
embodiment of portions of a trigger/marker module. 
20 Figure 4B is a schematic diagram illustrating generally one embodiment 

in which three electrode signals (e.g., RV, RM, LV) are used. 

Figure 5 is a timing diagram illustrating generally one embodiment of 
updating a heart rate variability marker. 

Figure 6 is a block diagram illustrating generally one conceptual 
25 embodiment of portions of an arrhythmia prediction module- 
Figure 7 is a block diagram illustrating generally one conceptual 
embodiment of portions of a preventive therapy control module. 

Figure 8 is a diagram illustrating generally one embodiment of a 
translation matrix. 
30 Detailed Description 

This invention relates to a cardiac rhythm management system that 
predicts when an arrhythmia will occur based upon the detection of conditioning 
events found to be statistically associated with the occurrences of arrhythmias. 
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A further embodiment uses the prediction to invoke a therapy to prevent or 
reduce the consequences of the arrhythmia. The present system may incorporate 
several features including: (1) distinguishing a degree of risk for occurrence of a 
future arrhythmia, (2) making predictions that apply to a well defined future 
5 time interval, (3) basing predictions on direct observations of arrhythmic 
triggers/markers (i.e., conditioning events) using conditional probabilities that 
these arrhythmias will occur given the past and present occurrences of the 
arrhythmias and the triggers/markers, (4) assessing a confidence in the accuracy 
of the prediction, (5) adapting its prediction and prevention capabilities to the 

1 0 individual patient, (6) assessing the expected effect of one or multiple preventive 
therapies in order to determine if the likelihood for future arrhythmias is 
increased or decreased by the proposed preventive therapies, (7) selecting a 
therapy based on this assessment of its expected effect, and (8) determining 
whether the selected preventive therapy combination should be invoked based on 

1 5 the magnitude, timing and confidence in the arrhythmia prediction and of the 
expected effect of the selected therapy. In one embodiment, the system also 
includes an external programmer device that can control, load and retrieve 
information from the implanted cardiac rhythm management device, and which 
can process and display information obtained from the implanted device. The 

20 present system provides techniques for predicting and preventing future cardiac 
arrhythmias, and such techniques can be used in combination with pacing, 
defibrillation, and other therapy techniques for treating presently existing cardiac 
arrhythmias. 

In the following detailed description, reference is made to the 
25 accompanying drawings which form a part hereof, and in which is shown by way 
of illustration specific embodiments in which the invention may be practiced. 
These embodiments are described in sufficient detail to enable those skilled in 
the art to practice the invention, and it is to be understood that the embodiments 
may be combined, or that other embodiments may be utilized and that structural, 
30 logical and electrical changes may be made without departing from the spirit and 
scope of the present invention. The following detailed description is, therefore, 
not to be taken in a limiting sense, and the scope of the present invention is 
defined by the appended claims and their equivalents. 

10 
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The present methods and apparatus will be described in applications 
involving implantable medical devices including, but not limited to, implantable 
cardiac rhythm management systems such as pacemakers, 
cardioverter/defibrillators, pacemaker/defibrillators, and biventricular or other 
5 multi-site coordination devices. However, it is understood that the present 
methods and apparatus may also be employed in unimplanted devices of the 
same sort, as well as in external monitors, programmers and recorders. 

Herein, the term "prediction" is used to mean a probability statement 
regarding whether or not an arrhythmia will occur at a future time. The term 
10 "trigger" refers to one or more cardiac events, such as depolarizations or other 
intrinsic heart activity signals, that may trigger an arrhythmia. Examples of such 
triggers include, by way of example, but not by way of limitation: sinus beats, 
premature sinus beats, beats following long sinus pauses; long-short beat 
sequences, R on T-wave beats, ectopic ventricular beats, and premature 
15 ventricular beats. The term "marker" refers to one or more quantities associated 
with at least one abnormal physiologic state (e.g., the occurrence of an 
arrhythmia) in which the quantity is obtainable from one or more 
electrophysiologic signals, or from one or more signals from one or more other 
sensors. Examples of such markers include, by way of example, but not by way 
20 of limitation: ST elevations, heart rate, increased or decreased heart rate, 

abnormal heart rate variability, late-potentials, or abnormal autonomic activity. 
Overview 

Figure 1 is a schematic drawing illustrating generally, by way of 
example, but not by way of limitation, one embodiment of portions of a cardiac 

25 rhythm management system 100 and an environment in which it is used. In 
Figure 1, system 100 includes an implantable cardiac rhythm management 
device 105, also referred to as an electronics unit, which is coupled by an 
intravascular endocardial lead 110, or other lead, to a heart 115 of patient 120. 
System 100 also includes an external programmer 125 providing wireless 

30 communication with device 105 using a telemetry device 130. In one 

embodiment, external programmer 125 includes a visual or other display for 
providing information to a user regarding operation of implanted device 105. 
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Catheter lead 110 includes a proximal end 135, which is coupled to device 105, 
and a distal end 140, which is coupled to one or more portions of heart 115. 

Figure 2 is a schematic drawing illustrating generally, by way of 
example, but not by way of limitation, one embodiment of device 105 coupled 
5 by leads 1 1 OA-C to heart 115, which includes a right atrium 200 A, a left atrium 
200B, a right ventricle 205A, a left ventricle 205B, and a coronary sinus 220 
extending from right atrium 200A. In this embodiment, atrial lead 11 OA 
includes electrodes (electrical contacts) disposed in, around, or near an atrium 
200 of heart 115, such as ring electrode 225 and tip electrode 230, for sensing 
1 0 signals and/or delivering pacing therapy to the atrium 200. Lead 1 10A 

optionally also includes additional electrodes, such as for delivering atrial and/or 
ventricular cardioversion/defibrillation and/or pacing therapy to heart 115. 

In Figure 2, a right ventricular lead HOB includes one or more 
electrodes, such as tip electrode 235 and ring electrode 240, for sensing signals 
15 and/or delivering pacing therapy. Lead HOB optionally also includes additional 
electrodes, such as coil electrodes 245A-B for delivering right atrial and/or right 
ventricular cardioversion/defibrillation and/or pacing therapy to heart 115. In 
one embodiment, system 100 also includes a left ventricular lead HOC, which 
provides one or more electrodes such as tip electrode 246 and ring electrode 247, 
20 for sensing signals and/or delivering pacing therapy. Lead HOC optionally also 
includes one or more additional electrodes, such as coil electrodes 248A-B for 
delivering left atrial and/or left ventricular cardioversion/defibrillation and/or 
pacing therapy to heart 115. 

In Figure 2, device 105 includes components that are enclosed in a 
25 hermetically-sealed enclosure, such as can 250. Additional electrodes may be 
located on the can 250, or on an insulating header 255, or on other portions of 
device 105, for providing unipolar pacing and/or defibrillation energy in 
conjunction with the electrodes disposed on or around heart 115. Other forms of 
electrodes include meshes and patches which may be applied to portions of heart 
30 115 or which may be implanted in other areas of the body to help direct electrical 
currents produced by device 105. For example, an electrode on header 255 may 
be used to stimulate local muscle to provide an alert/warning to the patient. In 
another example, an additional lead is used to provide electrodes associated with 

12 
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nerves or nerve ganglia, such as the vagus nerve, left or right stellate ganglion, 
carotid sinus nerve, or the fat pad over the atrioventricular node. The present 
method and apparatus will work in a variety of configurations and with a variety 
of electrical contacts or electrodes. 
5 Figure 3 is a schematic diagram illustrating generally, by way of 

example, but not by way of limitation, one embodiment of portions of device 
105, which is coupled to heart 115 and/or other portions of the patient's body. 
Figure 3 illustrates one conceptualization of various modules, which are 
implemented either in hardware or as one or more sequences of steps carried out 

10 on a microprocessor or other controller. Such modules are illustrated separately 
for conceptual clarity; it is understood that the various modules of Figure 3 need 
not be separately embodied, but may be combined and/or otherwise 
implemented, such as in software/firmware. 

In Figure 3, device 105 includes, among other things, a power source 

1 5 300, such as a battery. A communication module 305 includes a telemetry or 
other circuit by which implantable device 105 communicates with external 
programmer 125. A sensing module 310 senses intrinsic heart activity signals 
from one or more electrodes associated with heart 115. In one embodiment, 
sensing module 310 also senses other electrophysiological signals. A therapy 

20 module 315 provides therapy for treatment of present arrhythmias and 

prevention of future arrhythmias. In one embodiment, such therapy is provided 
at electrodes associated with heart 115 or portions of the nervous system such as, 
for example but without limitation, to sense other electrophysiological signals 
such as activity from sympathetic or parasympathetic members of the autonomic 

25 nervous system, or to sense blood temperature or blood flow. In various 

embodiments such therapy includes, among other things, pacing pulses, anti- 
tachyarrhythmia pacing (ATP), defibrillation shocks, etc. 

In Figure 3, sensing module 320 includes, among other things, one or 
more sensors, such as an accelerometer, acoustic sensor, respiration and/or stroke 

30 volume sensor (e.g., using transthoracic impedance), time and/or date of detected 
arrhythmia(s), cardiac displacement or blood vessel dimensions (e.g., using 
ultrasonic transit time measurements), or blood electrolyte levels (e.g., using ion- 
selective membranes). Sensing module 320 also includes interface circuits that 
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receive control signals and preprocess the sensor signal(s). Device 105 also 
includes a control circuit, such as a microprocessor or other controller 325, 
which communicates with various peripheral circuits via one or more nodes, 
such as bus 330. Controller 325 includes various functional blocks, one 
5 conceptualization of which is illustrated in Figure 3. 

In one embodiment, controller 325 includes a bradyarrhythmia control 
module 335 that detects bradyarrhythmias based at least in part on one or more 
signals obtained from sensing module 310. Bradyarrhythmia control module 
335 also provides one or more physiologically appropriate anti-bradyarrhythmia 

1 0 therapies to treat presently existing bradyarrhythmias. A tachyarrhythmia 
control module 340 detects tachyarrhythmias based at least in part on one or 
more signals obtained from sensing module 310. Tachyarrhythmia control 
module 340 also provides one or more physiologically appropriate anti- 
tachyarrhythmia therapies to treat presently existing tachyarrhythmias. 

1 5 For predicting and preventing arrhythmias, controller 325 also includes 

trigger/marker module 345, arrhythmia prediction module 350, and preventive 
therapy control module 355. Trigger/marker module 345 detects one or more 
triggers and/or markers based at least in part on signals received from sensing 
module 310 and/or sensing module 320. Arrhythmia prediction module 450 

20 predicts the likelihood of future arrhythmias using probability calculations based 
on trigger/marker information received from trigger/marker module 345. 
Preventive therapy control module 355 selects the most appropriate therapy (or 
combination of therapies) for preventing the future arrhythmia from a set of 
available preventive therapies. Preventive therapy control module 355 also 

25 triggers the delivery of such therapy after determining if the probability of 
arrhythmia, computed by arrhythmia prediction module 350, and the expected 
outcome of the selected therapy warrants administration of the therapy by 
therapy module 315. Thus, device 105 predicts and prevents future arrhythmias, 
as described more particularly below. In one embodiment, the prediction and 

30 prevention of future arrhythmias is deactivated if one or more present 

arrhythmias (e.g., ventricular tachyarrhythmia (VT) or ventricular fibrillation 
(VF)) are present and are being treated by techniques, known to one skilled in 
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the art, using tachyarrhythmia control module 340 and/or bradyarrhythmia 
control module 335. 

Example of Detecting Arrhythmogenic Trigger(s)/Marker(s) 

Figure 4A is a block diagram illustrating generally, by way of example, 
5 but not by way of limitation, one conceptual embodiment of portions of 
trigger/marker module 345. In this embodiment, trigger/marker module 345 
includes beat classification module 400, detection processing module 405, and a 
trigger/marker data bank such as trigger/marker list 410. Trigger/marker module 
345 recurrently examines signals from sensing module 310 and/or sensing 

10 module 320 and detects the presence, timing, and (if appropriate) magnitude of 
triggers/markers. These detections, timings, and magnitudes are output to 
detection processing module 405 either for each heartbeat, or corresponding to a 
time period encompassing multiple heartbeats. 

Beat classification module 400 examines signals from variously located 

1 5 electrodes. On a beat-to-beat basis, beat classification module 400 distinguishes 
between depolarization events that spread over the heart normally from those 
that spread out over the heart abnormally. Normal beats refer to those beats that 
are most prevalent and physiologically sound in the patient, even though the 
patient's most prevalent beat type might not be considered normal when 

20 compared to the propagation of a depolarization in a healthy heart. In one 
embodiment, sensing module 310 provides and beat classification module 
analyzes signals from: a right atrial sensing electrode (RA), a right ventricular 
sensing electrode (RV), a left ventricular (LV) sensing electrode, a right sided 
morphology (RM) signal between defibrillation electrodes in the right 

25 ventricular and superior vena cava and/or the metallic shell of the device, a left 
sided morphology (LM) signal between a left ventricular electrode and one or 
more of the other electrodes. Several example methods of operating beat 
classification module 400 to classify heart beats as normal or abnormal are 
described below by way of example, but not by way of limitation. 

30 A first method to distinguish normal and abnormal beats detects 

depolarizations on a plurality of electrode signals for each beat. In one example 
of this method, an RV depolarization indicates that a beat has occurred. The 
time intervals between detection of the same depolarization at the various 
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electrodes (RA, RV, LV, RM, LM) reflect the pattern and timing with which this 
depolarization propagates over the heart. Figure 4B is a schematic diagram 
illustrating generally, by way of example, but not by way of limitation, one 
embodiment in which 3 electrode signals (e.g., RV, RM, LV) are used. The 
5 sensed depolarization on each electrode signal resets 3 interval timers (total of 9 
timers) that measure time intervals initiated by the resetting event. The sensed 
depolarization on each electrode signal also latches the current value in 3 timers, 
each of which is reset by a different electrode signal, as illustrated in Figure 4B. 
In this way, each time a depolarization passes over the heart, the set of 9 time 

10 intervals is obtained. These time intervals represent the pattern with which the 
depolarization reached different electrodes. Similarly, if 4 signals are used, then 
16 intervals exist, and if 5 signals are used, then 25 intervals exist, etc. If the set 
of intervals for a particular beat fall outside of limits deemed normal, then that 
beat is classified as abnormal. In one embodiment, the abnormal limits vary 

1 5 with the heart rate. The abnormal limits may be programmed into the device as 
either population-based or patient-specific values. 

A second method to distinguish normal beats from abnormal beats 
compares the morphology of one or more available signals (e.g., RV and/or RM 
and/or LM signals) with corresponding template morphologies for a normal beat. 

20 In one embodiment, an RV depolarization detection provides a triggering event. 
For each RV depolarization, a correlation coefficient between the selected 
morphology signal and a corresponding template morphology is computed. If 
the coefficient is smaller than a predetermined threshold value, then the current 
beat is classified as abnormal. In one embodiment, different template 

25 morphologies are used for different ranges of heart rates. Also, the threshold 
values may be programmed into the device using either population-based or 
patient-specific values. 

A third method to distinguish normal and abnormal beats extracts 
features from one or more available signals and compares the extracted features 

30 with similar features for normal beats. Possible features include, by way of 
example, but not by way of limitation, QRS duration, R-wave amplitude, QT 
interval, and T-wave amplitude. For a particular beat, if the detected value of 
one or more features falls outside limits deemed normal, then the beat is 
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classified as abnormal. In one embodiment, the normal limits vary with heart 
rate. Also, the normal limits may be programmed into the device either as 
population-based or patient-specific values. 

In addition to these methods of classifying beats, other methods of 
5 identifying or classifying beats will also be suitable. Also, beat classification 
module 400 can also be used with lead configurations having a different number 
of electrodes or electrodes located in different locations. After beat classification 
module 400 classifies a detected heart beat, detection processing module 405 
performs the processing required for trigger/marker detection. 

1 0 In one embodiment, detection processing module 405 uses the beat-to- 

beat intervals and morphological data extracted from the signals by beat 
classification module 400 and also extracts any additional morphological 
measures required for trigger/marker detection. Detection processing module 
405 is also capable of examining activity, respiration, and stroke volume signals 

15 received from one or more acceleration, acoustic and/or impedance sensors and 
provided by sensing module 320. 

Device 105 includes a trigger/marker list 410, which lists the 
triggers/markers to be detected. List 410 includes members, for example, the 
types and natures of the particular trigger/marker sought to be detected. In one 

20 embodiment, each member of the list includes a corresponding predetermined 
comparison value. In one example, list 410 includes a member "long QRS 
duration" trigger/marker. For that member, list 410 also includes a comparison 
value indicating which QRS durations are long. In another example, list 410 
includes a particular "short-long-short sequence" of intervals between beats 

25 classified as normal. For that member, list 410 also includes corresponding 
comparison values establishing the criteria for long and short intervals in the 
sequence. Alternatively, list 410 includes members that include one or more 
beats classified as abnormal. 

In one embodiment, the type and nature of the triggers/markers in list 410 

30 are obtained from research in the relevant clinical population. Several different 
types of triggers/markers exist. For example, a first basic type of 
triggers/markers, is referred to as "sequence type" triggers/markers. Sequence 
type triggers/markers are based on sequences of one or more beats. Examples of 
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sequence type triggers/markers include, by way of example, but not by way of 
limitation, premature normal beats, premature abnormal beats, delayed normal 
beats, delayed abnormal beats, long-short intervals between beats, short-long 
intervals between beats, and multiple (couplets, triplets, etc) abnormal beats. 
5 Detection processing module 405 detects and counts the number of times each of 
the specified sequence type triggers/markers are present during a particular 
observation period. 

An example of a second basic type of triggers/markers is referred to as 
"value type" triggers/markers. Value type triggers/markers are based on metrics 

1 0 from the available signals, which are compared to predetermined comparison 
values to determine whether the trigger/marker is present during a particular 
observation period. Examples of value type triggers/markers based on 
morphology intrinsic heart activity signals include, by way of example, but not 
by way of limitation, QRS duration, ST magnitude, QT interval, and R-wave 

1 5 amplitude. In one embodiment, detection of morphology-based value type 

triggers/markers is based only on consideration of one or more beats classified as 
normal. In one example, a trigger/marker is based on the most recent beat. In 
another example, the trigger/marker is based on an average of morphology 
values over some number of previous normal beats. In a further example, the 

20 trigger/marker is based on one or more abnormal beats. Other examples of value 
type triggers/markers are not based on the heart's electrical activity. Such value 
type triggers/markers include, by way of example, but not by way of limitation, 
present respiratory rate, a position in the respiratory cycle at which the beat is 
detected, tidal volume of a respiratory inhalation or exhalation, cardiac stroke 

25 volume, present activity level of the patient, or the current position within the 
diurnal cycle. 

An example of a third basic type of triggers/markers is referred to as 
"history type" triggers/markers. In one embodiment, the presence or absence of 
history type triggers/markers is determined based on one or more possible 
30 computational analyses of historical data obtained from multiple beats or from 
multiple observation periods. Examples of history type triggers/markers include, 
by way of example, but not by way of limitation, a percentage of abnormal beats 
detected during an observation period, a percentage of premature or ectopic beats 
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detected during an observation period, heart rate variability during an 
observation period (e.g., 5 minutes), and the presence of altemans (i.e., cyclic 
variations over a plurality of cardiac cycles, e.g., T-wave alternans or QRS 
alternans) during an observation period. Other examples of history type 
5 triggers/markers include, by way of example, but not by way of limitation, short- 
term and long-term averages of the beat-to-beat data, trends in such averages 
(e.g., increasing or decreasing ST elevation, increasing density of abnormal or 
premature beats, etc.), and trends or periodicities in sensor values obtained from 
sensing module 320. 

10 In one embodiment, detection processing module 405 outputs a set of 

detection values, denoted D D 2 , D 3 , ... , D N , corresponding to the N members 
in list 410 and where each detection value, D t , corresponds to a single member of 
list 410. In one embodiment, a detection value Dj is set to 0 if the corresponding 
trigger/marker was not detected during the most recent observation period, and is 

1 5 set to 1 if the trigger/marker was detected during the most recent observation 
period. 

Timing of Arrhythmia Predictions and Trigger/Marker Detection 

Not all trigger/marker detections are updated on a beat-to-beat basis 
because some triggers/markers (e.g., history type triggers/markers) generally 

20 require information obtained from multiple beats. In one embodiment, the 
detection of triggers/markers is performed asynchronously from arrhythmia 
predictions. Several techniques can be used for this purpose. In one technique, 
each trigger/marker detection value corresponds to a predetermined observation 
period during which time its value (e.g., present, or not present) is determined. 

25 For many triggers/markers, the appropriate observation period is equal to a 

predetermined basic time period (BTP). For instance, if the basic time period is 
2 minutes, then device 105 tests for the occurrence of such triggers/markers in 
the last 2 minutes. However, detection observation periods are not necessarily 
the same for each trigger/marker. In one embodiment, for example, a heart rate 

30 variability trigger/marker is based on an observation period of approximately 5 
minutes. 

Figure 5 is a timing diagram illustrating generally, by way of example, 
but not by way of limitation, a heart rate variability trigger/marker, D 7 , which is 
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updated using a 5 minute observation period. In this example, D J = 1 during 
time periods 500 when heart rate variability is below its detection threshold and 
Dj = 0 during other time periods 505 when heart rate variability is above its 
detection threshold. Figure 5 also depicts a "long-short sequence" trigger, D 29 
5 which is updated at 2 minute intervals. In this example, D 2 = 1 during time 
periods 510 when the long-short sequence has been detected in the preceding 2 
minute observation period, and D 2 = 0 at other times. In this way, particular 
trigger/marker detection values are updated independently of other 
trigger/marker detection values in the set of trigger/marker detection values, D, 
1 0 as sufficient data for making the particular trigger/marker detection becomes 
available. 

In one embodiment, predictions of future arrhythmias are made at 
approximately regular time intervals 415, such as at the BTP. In one 
embodiment, the BTP is based on a fixed number of beats (e.g., 60 beats). In 

15 another embodiment, the BTP is based on the number of beats closest to a fixed 
duration (e.g., the number of beats just exceeding 1 minute, such that the BTP 
remains approximately constant. In a further embodiment, the time duration of 
the BTP is adjusted over a period of time during which device 105 is being used, 
such that predictions of future arrhythmias are made more or less frequently at 

20 times when the predicted probability for arrhythmias is higher or lower, 
respectively. 

Other Techniques for Timing Predictions of Arrhythmias 

According to one aspect of the present arrhythmia prediction and 
prevention techniques, the prediction provides a good match between the time 

25 period covered by the prediction and the requisite times for one or more potential 
preventive therapies. For example, a prediction of a 99.9% chance that the 
patient will have an arrhythmia within the next 100 years is highly accurate but 
worthless because the prediction does not relate to when the prevention therapy 
should be delivered or which prevention therapy should be provided. Similarly, 

30 a prediction of a 50% chance that the patient will have an arrhythmia within the 
next 2 years is of limited value because it relates only to long-term prevention 
therapies such as selecting patients for device implant or for long-term drug or 
pacing therapy. By contrast, a prediction of a 50% chance that the patient will 
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have an arrhythmia within the next 30 seconds provides great value in an 
arrhythmia prediction and prevention scheme in an implanted device provided 
that the decision to invoke a preventive therapy and the requisite time of action 
for that therapy can be carried out quickly enough to significantly reduce the 
5 chance of the future arrhythmia. In this example, the value of the prediction 
would be quite limited if the available preventive therapy required, for example, 
20 minutes to take effect. 

According to another aspect of the present arrhythmia prediction and 
prevention techniques, the prediction frequency adequately relates to the time 
10 period for which the prediction is valid. For example, it makes little sense for a 
prediction covering 1 month to be made once every 5 seconds or for a prediction 
covering 5 seconds to be made once every month. Thus, the present arrhythmia 
prediction and prevention techniques balance between the time period covered 
by a prediction, the frequency with which it is made, and the time required for 
the preventive therapy to take effect. 

In one embodiment, device 105 does not necessarily make the same set 
of predictions each time predictions are made. One method to do so for the 
device to make one set predictions for each BTP that cover times on the order of 
the BTP. Then at BTP multiples that correspond to longer intervals (e.g. 20 
minutes, 1 hour, 1 day, 1 week, 1 month), a corresponding sets of additional 
predictions are made which cover the correspondingly longer time intervals. The 
scheduling of the predictions, their corresponding time intervals are selected, 
according to the above argument, according to the times of action for the 
prevention therapies available to the device 105. 
Arrhythmia Prediction Example 

Figure 6 is a block diagram illustrating generally, by way of example, but 
not by way of limitation, one conceptual embodiment of portions of arrhythmia 
prediction module 350. In one embodiment, arrhythmia prediction module 350 
includes an arrhythmia probability calculation module 600, a conditional 
probability data bank such as a list of conditional probabilities 605, 
trigger/marker use data bank such as trigger/marker use list 610. In a further 
embodiment, arrhythmia prediction module 350 also includes an adaptive 
processing module 615. In one embodiment, conditional probability list 605 and 
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trigger/marker use list 610 each have members corresponding to the members of 
trigger/marker list 410, as explained below. 

In one embodiment, arrhythmia prediction module 350 includes an input 
that receives the detection values {D h . . . D N ) provided by trigger/marker 
5 module 345. In a further embodiment, arrhythmia prediction module 350 also 
includes input(s) that receive one or more "arrhythmia detected" signals 
provided by conventional arrhythmia detection modules of device 105 (e.g., 
signals from bradyarrhythmia control module 335 and tachyarrhythmia control 
module 340). These "arrhythmia detected" signals indicate, among other things, 

10 the presence or absence of one or more present bradyarrhythmias or 

tachyarrhythmias. Arrhythmia prediction module 350 outputs an arrhythmia 
prediction to preventive therapy control module 355, which, in turn, bases 
delivery of preventive therapy on the arrhythmia prediction. 

One aspect of the present arrhythmia prediction techniques treats the 

1 5 occurrence of an arrhythmia during a particular time period as a random event. 
In this embodiment, an arrhythmia prediction includes a probability calculation 
estimating the probability that an arrhythmia will occur within a specified period 
after the prediction. An example of such a prediction is a 50% probability that 
an arrhythmia will occur within the next 2 minutes. This method of prediction 

20 includes both a degree (e.g., 50%) and a well defined time period during which 
the prediction is applicable (e.g., 2 minutes), as opposed to merely indicating that 
the patient is currently at risk. 

In one embodiment, the arrhythmia prediction output from arrhythmia 
prediction module 350 includes a set of one or more arrhythmia probability 

25 assertions or statements. Each probability statement includes both a magnitude 
of the probability and a specified future time period associated therewith. In one 
embodiment, each probability statement also identifies which trigger, marker, or 
combination of triggers) and/or markers) contributed to its magnitude. In a 
further embodiment, the time period covered by each probability statement (i.e., 

30 the time period over which each probability statement is valid) is determined by, 
among other things, the scheduled prediction frequency (e.g., predictions made 
at 1 minute intervals cover a 1 minute period, etc.). 
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One conceptual approach to prediction recognizes that the rate at which 
arrhythmia events randomly occur changes relatively slowly compared to the 
time period over which predictions are made (i.e., the prediction frequency). 
Because future arrhythmia rates cannot be measured, predictions are based on 
5 one or more previous and/or present estimates of the arrhythmia rate assuming 
that the arrhythmia rate doesn't change significantly during the time period over 
which the prediction is made. The probability that a random arrhythmia event 
will occur within the specified prediction time period can be calculated based on 
the present arrhythmia rate. For simplicity, the occurrence of arrhythmia events 
1 0 can be understood as the result of a Poisson process with a rate, R. The 

probability, P, that one or more arrhythmias will occur within a next time period, 
T t is P = (1 - e '* T ) . The accuracy of P depends on the degree to which the 
estimate for R is correct, the rate R is truly constant during the time period T, and 
whether the underlying process is truly a Poisson random process. While such a 
1 5 conceptual approach may be helpful to understand the context of certain of the 
present prediction techniques, in one embodiment arrhythmia prediction module 
305 does not actually formulate an estimate of the underlying arrhythmia rate. 
Instead, the arrhythmia rate is reflected in the presence or absence of the detected 
triggers/markers, which are used by arrhythmia probability calculation module 
20 600 to formulate a probability for a future arrhythmia, as explained below. 

For example, consider the 7th trigger/marker of trigger/marker list 410. 
In one embodiment, the detection value D 7 is either one or zero depending on 
the respective presence or absence of this trigger/marker during the basic time 
period equal to T. Further, consider that the rate of arrhythmias Rj associated 
with this trigger/marker is zero when the trigger/marker is absent and that the 
rate R 7 has a nonzero value when the trigger/marker is present. The contribution 
of the 7th trigger/marker to the arrhythmia probability is D j{\ - o' RT \ or 
simply DjxCj where CP j is the conditional probability for the arrhythmia 
given that D } is present. The total probability for an arrhythmia, P, which 
includes the contributions from all trigger/markers, is computed as P = D } CP 
l +D 2 xCP 2 +D y xCP 3 +. . . -KD N xCP N , for the case where there is a detection value 
and a conditional probability for each member of the trigger/marker list 410. 
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In one embodiment, the conditional probabilities CP are obtained 
empirically from observations made in a patient population. In another 
embodiment, the conditional probabilities CP are seeded with population-based 
values, but are later adapted to the individual patient based on empirical 
5 observations of that patient. Because the conditional probabilities CP are 
empirically determined, the underlying process need not be exactly Poisson 
distributed, and the rate need not be completely unchanging during the time 
period T. The empirical estimations for the conditional probabilities CP 
incorporate deviations from these assumptions, although their predictive 

10 accuracy may be reduced. In one embodiment, a trigger/marker list 410 and a 
population-based conditional probability list 605 is loaded into device 105 
before, during, or after device 105 is implanted in the patient. 

In one embodiment, arrhythmia prediction module 350 also includes a 
trigger/marker use list 610, which provides a set of usage weights or flags, each 

15 usage weight or flag associated with a corresponding trigger/marker of 

trigger/marker list 410. According to one aspect of the present techniques, it 
may be desirable to exclude one or more of the trigger/marker detections from 
the total probability computation, such as, as the result of learning patient- 
specific traits by adaptive processing module 615. In trigger/marker use list 610, 

20 each usage weight or flag indicates how that trigger/marker detection value and 
conditional probability enters into the arrhythmia probability computation. In 
one embodiment, usage flags provide a binary indication of whether or not that 
particular trigger/marker detection value and conditional probability enters into 
the arrhythmia probability calculation. In another embodiment, usage weights 

25 provide a degree to which the detection value and conditional probability enters 
into the arrhythmia probability calculation. In such an embodiment, the 
arrhythmia probability calculation may be normalized, based on the values of the 
usage weights, such that the arrhythmia probability ranges between 0 and 1. 
Determining Triggers/Markers and Conditional Probabilities from a Population 

30 In one embodiment, a relevant patient population is used to obtain initial 

or actual estimates for conditional probabilities and/or to select particular 
triggers/markers for active use in the present prediction techniques. Patient- 
specific and adaptive techniques are described later. One approach for 
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determining population-based parameters uses long-term data (e.g., from sensing 
module 310 and/or sensing module 320) from a representative subset of the 
clinically relevant population. This long-term data obtained from a plurality of 
patients in the clinically relevant patient population is referred to as the 
5 population database. Data in the population database is divided into time periods 
equal to the BTP. After each such time period, trigger/marker detection 
processing is performed over the entire population database for all 
triggers/markers. The presence or absence of each trigger/marker (i.e., the set of 
detection values, D) is noted after each such time period. The presence or 

10 absence of one or more arrhythmias is also noted after each such time period. 
One population-based estimate for a particular conditional probability associated 
with the 7th trigger/marker (i.e. CP, ) is simply CP, = [(total number of time 
periods with arrhythmias that were preceded by a time period with /^present) ) 
(total number of time periods with D t present)]. A similar estimate is made for 

15 all members in trigger/marker list 410. This embodiment uses time periods 

equal to the same duration (i.e., the BTP) for all triggers/markers. Alternatively, 
the population database is divided into equal-length time periods that are 
different between individual triggers/markers. For example, for each 
trigger/marker, the population database may be divided up into equal length time 

20 periods that correspond approximately to the time period covered by that 
particular trigger/marker. 

Determining Predictive Capability of Particular Tri ggers/Markers 
Not all possible triggers/markers may have predictive power, either in the 
clinically relevant patient population, or in the particular patient toward which 

25 the disclosed techniques are applied. The present techniques include a method of 
selecting those particular triggers/markers, in the trigger/marker list 410, that 
have useful predictive capability. One such method uses a Chi squared test 
based on the population database, but other statistical test methods may also be 
used. For example, in the Chi squared test, the presence or absence of 

30 arrhythmias and the available triggers/markers in each time period in the 

population database yields the following values: #BTP = total number of BTP's 
examined, #Z) + = total number of BTP's in which the trigger/marker was 
detected, #D" = total number of BTP's in which the trigger/marker was not 
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detected, #A + = total number of BTP's in which arrhythmia was detected, #A~ = 
total number of BTP's in which arrhythmia was not detected, #D*A + = number 
of BTP with the trigger/marker detected followed by a BTP with arrhythmia, # 
D + A' = number of BTP with the trigger/marker detected followed by a BTP 
5 without arrhythmia, # DA* = number of BTP without the trigger/marker detected 
followed by BTP with arrhythmia, #D'A' = number of BTPs without the 
trigger/marker detected followed by BTP without arrhythmia. 

A statistical test determines whether the observed occurrences of #D + A* y 
#D*A\ UD'A*, and WA' are different from those that would be expected if the 
1 0 trigger/marker did not have predictive power. One example such test computes 
the following sum: 

Sum =[ M+A+-#A + x UDl I #RTP] 2 + [ Mix un + / urtp ] 2 

#A + x #ZT ) #BTP M x W + I #BTP 

15 

+ [ #ftVr-#/Tx m /#BTP] 2 + [ # DA' - #/4' x UP' j MRTP ] 2 
#A + x m I UBTP M x m I #BTP 

This represents the Chi Square value (1 degree of freedom) that tests if 
20 the arrhythmia is associated with the trigger/marker. If the sum exceeds 3.84, for 
example, we are 95% confident that the trigger/marker provides predictive 
capability for the arrhythmia. A more complete multi-variate statistical analysis 
would provide stronger conclusions about the trigger/marker's role when 
considered together with other triggers/markers having independent predictive 
25 power. Strictly speaking, the form for the above probability computation is 
based on an approximation, i.e., that the predictive capability of the individual 
triggers/markers are independent from each other. 
Patient Specific Adaptive P rocessing Example 

In one embodiment, device 105 includes a patient specific adaptive 
30 processing module 615, as illustrated in Figure 6. In one example, adaptive 
processing module 615 modifies conditional probability list 605 and 
trigger/marker use list 610. The present techniques recognize that it is unlikely 
for some set of one or more triggers/markers to be so powerful at predicting 
arrhythmias that they always worked in all patients. The predictive capability of 
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a particular set of one or more triggers/markers will most likely vary between 
different patients. Population-based average values may be suboptimal for a 
particular patient. In one embodiment, device 105 examines long-term data 
specific to the individual patient in which it is implanted. Thus, one aspect of 
5 the present techniques allow device 105 to adapt to an individual patient to 
improve the accuracy and confidence in arrhythmia predictions over a time 
during which device 105 is being used. 

One embodiment extracts patient-specific parameters (e.g., conditional 
probabilities and/or weights) analogously to the techniques described above that 

1 0 extract population-based parameters. In this embodiment, the members of the 
trigger/marker list 410 are initially established based on clinical research in a 
population of interest, and are loaded into device 105. After implantation of 
device 105 in a particular patient, as ongoing data is acquired from the patient, 
adaptive processing module 615 manages trigger/marker use list 610 by 

1 5 determining which members of trigger/marker list 410 provide predictive 
capability for arrhythmia prediction in the particular patient. 

One such technique initially sets the members of the trigger/marker use 
list 610 to "do not use." At some point after implant, adaptive processing 
module 615 examines each BTP to detect the presence or absence of (1) 

20 arrhythmias and (2) the triggers/markers in trigger/marker list 410 to form a 
patient-specific database. Device 105 recurrently performs the above Chi 
squared test to determine (to within a predetermined level of confidence) if each 
member of trigger/marker list 410 provides predictive capability in the particular 
patient. If a member does provide such predictive capability, the usage value for 

25 that member is set to "use" in trigger/marker use list 610. Such a technique only 
allows an arrhythmia prediction based on a particular trigger/marker after 
sufficient confidence exists that its predictive capability in the particular patient 
is sufficiently high. 

Another such technique initially sets the members of the trigger/marker 

30 use list 610 to "use," such as where population based data indicates the 

likelihood that these members of the trigger/marker use list are likely to have 
sufficient predictive capability. In this embodiment, the corresponding usage 
values are set to "do not use" if the patient-specific database demonstrates that 
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the trigger/marker does not provide adequate predictive capability in the 
particular patient despite its good performance in the population. Such a 
technique stops using the trigger/marker to make predictions after sufficient 
confidence exists that its predictive capability in the particular patient is 
5 inadequate. 

For patient-specific adaptive processing, the ratios ( R* p = #D + A + ) # 

Z) + ) and (R~ P = UD 'A* ) #D~) from the population database (hence the P 

subscripts) are used to compute the following sum: 

Sum =[ # D+ A * - ffVx # D+] 2 + [ #D + A ' - (1 - ft*r) x# /T] 2 
10 /T/> x#£> + (1-/? + p)@ #D + 

+ [ O IL A l - fl>x# fl'] 2 + [ #DA - - (! - /?» x #D '] 2 
R'px #D- (UR'p)x #D- 

where the values of #ZT, #D*A*, etc. are obtained from the patient specific 
1 5 database. This sum, which is a Chi squared value, is compared to another 
predetermined threshold value. If the sum exceeds the threshold, then that 
trigger/marker does not predict the arrhythmia in the specific patient in the same 
manner as it does in the population. This could occur either if the trigger/marker 
provides lesser predictive power than in the population, or also if the 
20 trigger/marker provides greater predictive power than in the population. 

Consequently, the corresponding usage flag in trigger/marker use list 610 is only 
set to "do not use" if the patient-specific test also demonstrates that the 
trigger/marker does not provide predictive power in the particular patient. In one 
embodiment, these threshold values (population-based and/or patient-specific) 
25 are programmable. In a further embodiment, these threshold values (population- 
based or patient-specific) are programmed to either the same value, or to 
different values. 

The predictive capability of a particular trigger/marker could change over 
time in a particular patient. In one embodiment, adaptive processing module 615 
30 recurrently sets and/or resets the usage values in the trigger/marker use list 610 
over time while device 105 is being used. In this embodiment, testing for 
patient-specific predictive capability of each trigger/marker is carried out 
recurrently. If a particular trigger/marker is not presently in use, then the test 
result (i.e., the sum described above) would be required to exceed a first 
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threshold value to change the usage value of the trigger/marker to <4 use." If the 
trigger/marker is presently in use, then the test result (i.e., the sum described 
above) would be required to exceed a second threshold value to change the usage 
value of the trigger/marker to "do not use." These threshold values can either be 
5 the same, or alternatively can be different to provide hysteresis that would 
reduce the number of spurious transitions. 

In one embodiment, adaptive processing module 615 also provides and 
adaptively modifies conditional probabilities in conditional probability list 605 
based on observations in the specific patient. For a particular patient, a 
10 trigger/marker's conditional probability may vary from the population-based 
value. Conditional probability list 605 is initially seeded with population-based 
conditional probabilities. However, as device 105 acquires data from the 
particular patient, patient-specific conditional probabilities are used, as described 
below. 

1 5 In one embodiment, adaptive processing module 615 implements the 

patient-specific conditional probabilities by recurrently estimating the 
conditional probability for each trigger/marker analogously to their population- 
based estimation described above. However, patient-specific data is used rather 
than population-based data. For each member of the trigger/marker list 410, 

20 device 105 estimates the patient-specific conditional probability as CP Iest = # 
Dj + A + ) # where # Z>/ is number of BTP's in which the 7th trigger/marker 
was detected, and # D^A* is the number of BTP in which (1) trigger/marker I 
was detected, and (2) the BTP was followed by a BTP with an arrhythmia. 

Such recurrent reforecasting of conditional probabilities also includes a 

25 second step of determining when the patient-specific values should be used 

rather than the population-based values. Counters storing #D^A + and #£>/ are 
initially reset to zero. Whenever a BTP containing trigger/marker l is detected, 
device 105 formulates a new estimate for the conditional probability CP Iest by 
incrementing #Z)/ and #Dj + A + if the following BTP includes an arrhythmia. 

30 If the initial value for CP, is also correct for the particular patient, then 

after observing # £>/ BTP's with trigger/marker 7, one would expect (# £)/ x 
CP j ) cases where an arrhythmia was found. More or fewer arrhythmias may 
result from chance alone. The standard deviation for the number of expected 
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arrhythmias is { #£/ x CP, (1- CP,) } ,/2 . A 95% confidence interval would 
include about 1.96 times this standard deviation above and below the expected 
values. Thus, if the initially entered values for CP, were also valid for this 
particular patient, then there exists a 95% confidence that the number of 
5 arrhythmias observed after # Z>, + occurrences of trigger/marker / is between [# 
D; x CP, - 1.96 {#d; x CPjX (1- CP,) } m ] and [#D/ x CP, + 1.96 { #£/ x 
CP, x (1- CP,) } 1/2 ]. Expressed as percentages of #Z)/, this yields the following 
confidence interval for the conditional probability: [CP, - 1 .96 x {CP, x (1- CP,) 
/#Z), + } ,/2 ] and[CP, +1.96x { CP,x (1- CP,) / #Z), + } 1/2 ] 

10 If the estimate for the conditional probability CP Jest falls outside this 

range, there exists a 95% confidence that CP Jest is a better estimate for this 
particular patient than CP, . In that case, the estimated value for the conditional 
probability is used for subsequent predictions of future arrhythmias rather than 
the population-based values. Because the patient's conditional probabilities may 

1 5 change over time, the same procedure can be repeated recurrently by using the 
new conditional probability as the given value and forming another new estimate 
based on a new set of observations. In this way, the conditional probabilities are 
recurrently updated to more appropriate values as device 105 is being used. 

One alternate method of updating conditional probabilities recurrently 

20 updates the patient-specific conditional probability for each trigger/marker in 
such a way that the values eventually approach or track the true values for the 
specific patient. In this method, the conditional probabilities are seeded with 
population-based values. Then, each time a BTP containing one of the 
triggers/markers is observed, device 105 notes the presence or absence of an 

25 arrhythmia during the next BTP. If an arrhythmia is present, then the conditional 
probability for that trigger/marker is increased by a small step (e.g., CP, = CP, 
+ STEP * (1-CP,)). If the arrhythmia is not present, then the conditional 
probability for that trigger/marker is decreased by a small step, as illustrated 
below (e.g., CP, = CP, - STEP * (CP,)). 

30 If the patient's true conditional probability equals CP, , then the ratio of 

BTP's with arrhythmias absent to those with arrhythmias present should equal ( 
1 -CP, ) / CP,. For a conditional probability of 20%, for example, an average of 
4 BTPs without arrhythmias are expected for each BTP with an arrhythmia. In 
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this case, the conditional probability is decreased 4 times by a small amount (0.2 
x STEP) but increased 1 time by a large amount (0.8 x STEP), yielding no net 
drift in the CP f over time. If the patient's actual conditional probability was 
significantly higher than CP h then there would be a higher proportion of BTP' s 
5 with an arrhythmia. As a result, CP f would tend to increase until CP, approached 
the correct value. Similarly, if the actual value was too low, CP 7 would decrease 
until it approached the correct value. In one embodiment, the selection of STEP 
is small (~ 0.05) such that these changes are stable over time. 
Examples of Alternate Methods for Arrhythmia Prediction 

10 In one alternate embodiment, arrhythmia prediction module 350 

incorporates into an arrhythmia prediction the number of times a trigger/marker 
occurs during a BTP or other time period. An arrhythmia may occur each time a 
particular trigger/marker occurs. Repeated appearances of a trigger/marker may 
suggest a stronger relative predictive value. Consequently, multiple occurrences 

15 of a trigger/marker provide a higher arrhythmia probability. In this embodiment, 
each member in the conditional probability list 605 includes a set of one or more 
conditional probability values. In one example, the first value in the set reflects 
the conditional probability when the corresponding trigger/marker occurred once 
during the basic time period, the second set reflects the conditional probability 

20 when the corresponding trigger/marker occurred twice during the basic time 
period, etc. In this embodiment, the set of conditional probabilities may 
alternatively reflect ranges of the number of trigger/marker detections (e.g., 
between 3 and 5 detections, fewer that 8 detections, more than 3 detections, etc.). 
The population-based and patient-specific values for these sets of conditional 

25 probabilities are also obtained by extension of the techniques described above. 
The patient-specific adaptive processing for these sets of conditional 
probabilities is also obtained by extension of the techniques described above. 

This alternate embodiment also provides a different technique of 
trigger/marker detection. Instead of setting the detection value, Z> 7 , to either zero 

30 or one, the trigger/marker detection processing module 405 sets the detection 
value D f to the number of times the 7th trigger/marker occurred during the BTP. 
The arrhythmia probability calculation is computed as P = D f x P IDl +D 2 x 
CP 2D2 + D 3 x CP 3D3 + ... + D N x CP NJ)N , where CP 1K is the conditional 
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probability when trigger/marker / is present K times during the BTP, and N is the 
number of triggers/markers. 

A second alternate embodiment recognizes that if two (or more) different 
triggers/markers occur during the same BTP, the probability for arrhythmias may 
5 be different from the linear sum of the conditional probabilities when the 
trigger/marker occurred alone (i.e., without the other one or more 
triggers/markers). In this embodiment, detection processing module 405 outputs 
binary detection values for the members in the trigger/marker list 410, However, 
in this embodiment, arrhythmia prediction module 350 forms a number from 0 ( 

10 none of triggers/markers in trigger/marker list 410 detected) to 2 N -1 (all of the 
triggers/markers in trigger/marker list 410 detected) indicating which of the 
triggers/markers in trigger/marker list 410 were detected during the BTP, where 
N is the number of members in trigger/marker list 410. In this embodiment, the 
members of conditional probability list 605 and trigger/marker use list 610 do 

1 5 not correspond to the members of trigger/marker list 41 0. Instead, the members 
of conditional probability list 605 and trigger/marker use list 610 correspond to 
the various combinations of the 2 N -1 possible combinations of triggers/markers 
that could be detected. The arrhythmia probability calculation is computed as P 
= CP K , where A' is a number (from 0 to 2"-l) that reflects the current 

20 combination of detected triggers/markers (i.e., those triggers/markers having a 
detection value of 1). 

Automated Preventive Therapy Selection Processing Example 

Preventive therapy control module 355 automatically decides whether to 

invoke preventive therapy based on arrhythmia probability statements obtained 
25 from arrhythmia prediction module 350. Preventive therapy control module 355 

provides output signals to therapy module 315 to deliver pacing and/or shock 

and/or other arrhythmia prevention therapies. 

Figure 7 is a block diagram illustrating generally, by way of example, but 

not by way of limitation, one conceptual embodiment of portions of preventive 
30 therapy control module 355. In this embodiment, prediction scheduler 700 

schedules predictions of future arrhythmias. Preventive therapy decision module 

705 decides whether arrhythmia prevention therapy is warranted. Preventive 

therapy selection module 710 selects one or more appropriate prevention 
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therapies. Prevention activation module 715 activates the selected arrhythmia 
prevention therapy. Preventive therapy control module 355 also includes a 
prevention therapy list 720, and a trigger/marker vs. preventive therapy 
translation matrix 725 that relates the prevention therapies of preventive therapy 
5 list 720 to the triggers/markers used by arrhythmia prediction module 350 in 
predicting future arrhythmias. The various submodules in therapy control 
module 355 are illustrated as such for conceptual purposes only; alternatively, 
these submodules may be understood as being incorporated in arrhythmia 
prediction module 350 or elsewhere. 

1 0 In one embodiment, prevention therapy list 720 includes all the possible 

arrhythmia prevention therapies that device 105 can deliver to the patient. List 
720 is programmed into device 105 either in hardware, firmware, or software. 
The members of list 720 may be selected based, by way of example, but not by 
way of limitation on: the patient's lead and electrode configuration, the patient's 

15 dependence on pacing, functionality of the patient's AV node, etc. In one 

embodiment, each member of preventive therapy list 720 includes information 
fields (e.g., pacing rate, interval sequences, etc.) that establish characteristics of 
that particular arrhythmia prevention therapy. In one embodiment, preventive 
therapy list 720 includes immediate, short-term, intermediate-term, and/or long 

20 term preventive therapies. 

Immediate preventive therapies include, by way of example, but not by 
way of limitation: overdrive atrial pacing, such as for patients with functional 
AV nodes; demand atrial pacing, such as for patients with functional AV nodes; 
overdrive ventricular pacing; demand ventricular pacing; simultaneous or 

25 sequenced overdrive atrial and ventricular pacing; simultaneous or sequenced 
demand right ventricular and left ventricular overdrive pacing, such as for 
patients with a left ventricular pacing electrode; simultaneous or sequential right 
ventricular and left ventricular demand pacing, such as for patients with a left 
ventricular pacing electrode. 

30 Short-term preventive therapies include, by way of example, but not by 

way of limitation: repetitive overdrive atrial pacing, such as for patients with 
functional AV nodes; repetitive simultaneous or sequenced demand atrial and 
ventricular pacing; subthreshold cardiac stimulation, such as using defibrillation 
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lead electrodes; neural stimulation of the autonomic nervous system, such as for 
patients with autonomic stimulation leads; global pacing pulses implemented by 
providing low energy shocks via defibrillation electrodes; delivery of certain 
drugs. 

5 Intermediate-term preventive therapies include, by way of example, but 

not by way of limitation, an alert/warning for the patient and/or physician such 
as by providing an audible tone or other signal; delivery of certain drugs. 

Long-term preventive therapies include, by way of example, but not by 
way of limitation, diagnostic warnings for the physician, such as by transmitting 

10 diagnostic information from device 105 to external programmer 125 using a 
telemetry or other communication link; delivery of certain drugs. 

According to one aspect of the present system, each member of 
preventive therapy list 720 is associated with a required time of action, which 
includes one or more of a time for the therapy to become effective and/or a time 

1 5 after which the therapy is no longer effective. Accordingly, in one embodiment, 
the prediction scheduler 700 considers only those members of the preventive 
therapy list that can be expected to be effective within a time frame 
commensurate with the prediction time period. In other embodiments, 
information relating to the time frame in which a particular therapy is expected 

20 to be effective is incorporated into the translation matrix 725, and the prediction 
time period is added to the trigger/marker list 410. Other physiological variables 
can also be added to the list 410 to form a kind of physiological state vector that 
can be mapped to a particular member of the preventive therapy list 720, with the 
translation matrix incorporating information relating to the appropriateness each 

25 member of the therapy list for a given value of a physiological variable. 

In one embodiment, only one member of the preventive therapy list 720 
is invoked at any particular time. Combinations of different preventive therapies 
are also provided, but each such combination is treated as a separate entry in 
preventive therapy list 720. For example, "initiate overdrive ventricular pacing" 

30 is one member in preventive therapy list 720 while "initiate an alert/warning 
stimulus" is a different member in preventive therapy list 720. A combined 
therapy using both "initiate overdrive ventricular pacing" and "initiate an 
alert/warning stimulus" is treated as yet another entry in preventive therapy list 
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720. This simplifies the tasks of preventive therapy selection module 710 and 
prevention activation module 715. 

Preventive therapy selection module 710 selects an arrhythmia 
prevention therapy based on outputs from preventive therapy decision module 
5 705. If preventive therapy decision module 705 determines that the degree and 
confidence in the arrhythmia prediction warrant some preventive therapy, as 
discussed above, then preventive therapy selection module 710 selects a member 
of the preventive therapy list 720 to be invoked. In one embodiment, the 
selection of a prevention therapy is based on the set of trigger/marker detection 

10 values, D, upon which preventive therapy decision module 705 based the 
decision to provide arrhythmia prevention therapy. Translation matrix 725 
translates between the trigger/marker detection values D and selection of the 
appropriate arrhythmia prevention therapy from preventive therapy list 720. 

Figure 8 is a diagram illustrating generally, by way of example, but not 

1 5 by way of limitation, one embodiment of a translation matrix 725. In this 
embodiment, rows of translation matrix 725 correspond to the members of 
trigger/marker list 410. Columns of translation matrix 725 correspond to the 
members of preventive therapy list 720. In one embodiment, the cells in 
translation matrix 725 are set to either "Use," "Do Not Use," or "Don't Care." A 

20 "Use" entry indicates that the corresponding preventive therapy is expected to 
provide a beneficial effect to reduce the risk of the arrhythmia when the 
associated trigger/marker is present. A "Do Not Use" entry indicates that the 
corresponding preventive therapy is expected to result in a detrimental effect that 
may increase the risk of arrhythmia when the associated trigger/marker is 

25 present. The 'Don't Care " entry indicates that the corresponding preventive 
therapy is expected to neither increase or decrease the risk of arrhythmia when 
the associated trigger/marker is present. 

For the particular example illustrated in Figure 8, when the first member 
of trigger/marker list 410 is present, preventive therapies #1 and #6 are likely to 

30 be useful, but preventive therapy #4 should not be used because it is likely to be 
detrimental. The other possible preventive therapies are deemed to have no 
effect when trigger/marker #1 is present. 



35 



WO 01/24876 PCT/US00/26679 

In one embodiment, the values in translation matrix 725 are population- 
based values that are initially loaded into device 105. These population-based 
values are obtained for the population of interest based on observations of: the 
occurrences of arrhythmias; the presence of preceding triggers/markers; and the 
5 estimated or empirically derived effect of the particular preventive therapies on 
decreasing the probability of arrhythmia. 

In one embodiment, when preventive therapy decision module 705 
indicates that a preventive therapy is needed, there will be some set of detected 
triggers/markers that initiated the decision to invoke arrhythmia prevention 
10 therapy. Therapy selection module 710 considers those rows in translation 
matrix 725 that correspond to those particular triggers/markers that initiated 
arrhythmia prevention therapy. Selection of the appropriate preventive therapy 
is based on a score determined for each of the preventive therapies in preventive 
therapy list 720. The particular preventive therapy with the highest score is 
1 5 selected by preventive therapy selection module 710. 

In one embodiment, each member of preventive therapy list 720 is 
considered. For those rows that correspond to the presently detected 
triggers/markers, the number of "Use" entries is added to obtain the score. If any 
"Do Not Use" entry is present, however, the score is set to zero. In this way, 
20 each of the preventive therapies has a score between zero and the number of 
detected triggers/markers having a corresponding "Use" entry. Therapy 
selection module 710 selects that preventive therapy with the highest score. 
Where there is a tie, the preventive therapy appearing earliest in preventive 
therapy list 720 is used. Where all possible preventive therapies have a score of 
zero, then none of the available preventive therapies are deemed appropriate, and 
no preventive therapy is invoked. In one alternate embodiment, the score is also 
based on the conditional probability for each detected trigger/marker. In another 
alternate embodiment, a threshold value is included, such that the highest scoring 
preventive therapy must have a score that reaches (or exceeds) the threshold 
value before any preventive therapy is selected. 

Prevention activation module 715 uses the information contained in the 
preventive therapy list 720 and/or translation matrix 725 for activating delivery 
of the selected therapy by therapy module 315. Device 105 delivers the 
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activated therapy for a predetermined amount of time, which is also included in 
the preventive therapy list 720. The arrhythmia prevention therapy is delivered, 
by way of example, but not by way of limitation, by controlling the pacing 
and/or shock circuits or by storing information for later review (via telemetry or 
5 other communication to external programmer 125) by the physician or other 
user. 

Example of Scheduling Predictions 

Prediction scheduler 700 schedules the frequency with which arrhythmia 
predictions are made. As discussed earlier, the BTP used to detect a 

1 0 trigger/marker for prediction may differ from the time period covered by a 
particular prediction which, in turn, may differ from the time required for the 
preventive therapy to have an effect. Prediction scheduler 700 should schedule 
predictions in a rational manner, for allowing efficient operation. 

As described above, the arrhythmia probability statements can be 

1 5 conceptualized as having the following basic form: "In the next Y time period 
there is an X percent chance of arrhythmia due to the combination of 
triggers/markers D" In one embodiment, predictions are made at a maximum of 
2 minute intervals, and the BTP for detecting triggers/markers is also 2 minutes. 
In an alternate embodiment, however, the prediction frequency may be different 

20 from the BTP frequency. Intervals between predictions may be longer or shorter 
than the BTP, and may be variable or fixed. In one embodiment, the time period 
covered by the predictions (i.e., the values for Y) are either 2 minutes (i.e., 1 
BTP), 30 minutes (i.e., 15 BTPs), 1 day, 1 week, or 1 month. Predictions 
covering a 2 minute period are made at 2 minute intervals. Predictions covering 

25 a 30 minute period are made at 10 minute intervals. Predictions covering a 1 day 
period are made at 6 hour intervals. Predictions covering a 1 week period are 
made daily. Predictions covering a 1 month period are made weekly. Thus, in 
this embodiment, the predictions covering shorter time periods are made at a 
frequency corresponding to the time period covered by the prediction. 

30 Predictions covering longer time periods are made more often than the time 
period covered by the prediction. The timing of the predictions relative to the 
time covered by the prediction could be made in many other ways without 
departing from the techniques disclosed in this document. 
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Preventive Th erapy Decision 

Preventive therapy decision module 705 decides whether preventive 
therapy is warranted. In one embodiment, this task is performed at each 2 
minute interval (i.e., every BTP), at which time it considers all arrhythmia 
5 predictions existing at that time. In one embodiment, the total probability for 
arrhythmia for each of the prediction time periods is compared to corresponding 
threshold values. 

If the arrhythmia probabilities during all of covered time periods fail to 
exceed the corresponding threshold values, then no significant arrhythmia risk is 
10 deemed to exist and no preventive therapy is provided. If the arrhythmia 
probability during any of the covered time periods exceed the corresponding 
threshold value, then a significant arrhythmia risk is deemed to exist, and the 
arrhythmia prevention therapy is selected and activated as described above. 
Conclusion 

1 5 Although the above description described particular embodiments using 

an implanted cardiac rhythm management device including defibrillation 
capability, the techniques can also be used in other cardiac rhythm management 
systems including, without limitation, implanted or external pacemakers, or 
other acute or chronic cardiac care or monitoring devices. Moreover, although 

20 the techniques were described using those electrodes and sensors available in 
implanted cardiac rhythm management devices, different and/or additional 
sensing and/or stimulation electrodes may be used (e.g., for sensing or 
stimulating sympathetic or parasympathetic nerves or ganglion, or for sensing 
pH, p0 2 or K+ concentration in blood, etc.). 

25 It is to be understood that the above description is intended to be 

illustrative, and not restrictive. Many other embodiments will be apparent to 
those of skill in the art upon reviewing the above description. The scope of the 
invention should, therefore, be determined with reference to the appended 
claims, along with the full scope of equivalents to which such claims are entitled. 
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What is claimed is: 

1 . A method, comprising: 

detecting a conditioning event statistically associated with the occurrence 
5 of an arrhythmia in a patient's heart; and, 

predicting the occurrence of an arrhythmia within a specified prediction 
time period if an estimated arrhythmia probability exceeds a specified threshold 
value, wherein the estimated arrhythmia probability is computed from a 
conditional arrhythmia probability associated with the conditioning event that is 
1 0 derived from past observations of instances in which the conditioning event 
occurs alone or together with an arrhythmia within a specified time period. 



2. The method of claim 1 wherein the conditioning event is a marker event 
corresponding to detection of a physiological state statistically associated with 

1 5 occurrence of cardiac arrythmias. 

3. The method of claim 1 wherein the conditioning event is a trigger event 
that increases the risk of an arrhythmia occurring via a depolarization that serves 
as a source for the arrhythmia. 

20 

4. The method of claim 1 wherein the conditioning event is a specific 
morphology of a waveform representing electrical activity of the heart. 

5. The method of claim 1 wherein the conditioning event is a specific 

25 pattern of activation times of different areas of the heart as sensed by a plurality 
of electrodes. 



6. The method of claim 1 wherein the conditioning event is a specific 
sequence pattern of heartbeats with respect to time. 

30 

7. The method of claim 1 wherein the conditioning event is a value of a 
measured physiological variable. 
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8. The method of claim 1 wherein the conditioning event is a statistic based 
upon a history of occurrences of conditioning events. 

9. The method of claim 1 wherein the conditional arrhythmia probability 

5 CP is a ratio of the number of observed instances in which the conditioning event 
is followed by an arrhythmia within a specified basic time period to the total 
number of observed instances of the conditioning event. 

1 0. The method of claim 9 further comprising estimating a rate C at which 
10 the conditioning event occurs, and further wherein the estimated arrhythmia 

probability is calculated by the expression: 

estimated arrhythmia probability = CP x ( 1- e" 07 ). 

1 1 . The method of claim 1 wherein the conditional arrhythmia probability is 
15 calculated by the expression: 

l-e RT 

where T is a measure of the specified prediction time period, and R is an 
estimate of the rate at which arrhythmias occur while the conditioning event is 
present. 

20 

12. The method of claim 1 1 wherein the rate R is a ratio of: 1) the number of 
instances in which the conditioning event is followed by an arrhythmia within a 
specified basic time period, to 2) the length of the basic time period multiplied 
by the total number of basic time periods in which the conditioning event is 

25 observed. 

13. The method of claim 12 further comprising estimating a rate C at which 
the conditioning event occurs, and further wherein the estimated arrhythmia 
probability is calculated by the expression: 

30 (l-e RT )(l-e^ T ). 



14. The method of claim 1 further comprising assigning a culpability to a 
conditioning event for causing an observed instance of an arrhythmia when a 
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specified criterion is met, and further wherein the conditional arrhythmia 
probability CP is a ratio of the number of observed instances in which the 
conditioning event was deemed culpable for causing an arrhythmia, to the total 
number of observed occurrences of the conditioning event. 

5 

15. The method of claim 14 further comprising estimating a rate C at which 
the conditioning event occurs, and further wherein the estimated arrhythmia 
probability is calculated by the expression: 

estimated arrhythmia probability = CP x ( 1- e" CT ). 

10 

1 6. The method of claim 1 4 wherein the culpability criterion is the proximity 
in time between the conditioning event and the onset of an arrhythmia. 



17. The method of claim 14 wherein the culpability criterion is the 

1 5 magnitude of the detected conditioning event prior to the onset of an arrhythmia. 

18. The method of claim 14 wherein the culpability criterion is the frequency 
of occurrence of the conditioning event within a specific time period prior to the 
onset of an arrhythmia. 

20 

19. The method of claim 1 wherein a plurality of conditioning events 
statistically associated with the occurrence of an arrhythmia are detected and 
further wherein a composite estimated arrhythmia probability is compared with 
the threshold value in order to predict the occurrence of an arrhythmia, the 

25 composite arrhythmia probability being a combination of the estimated 
arrhythmia probabilities associated with each detected conditioning event. 

20. The method of claim 19 wherein the composite arrhythmia probability is 
calculated by adding the estimated arrhythmia probabilities associated with each 

30 detected conditioning event. 

2 1 . The method of claim 20 wherein the estimated arrhythmia probabilities 
associated with each detected conditioning event are calculated by multiplying 
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the conditional arrhythmia probability derived for a conditioning event by the 
number of detected occurrences of the event within a specified time period. 

22. The method of claim 1 9 wherein each conditional arrhythmia probability 
5 is a ratio of the number of instances in which the conditioning event is observed 

to be followed by an arrhythmia within a specified basic time period to the total 
number of observed instances of the conditioning event. 

23. The method of claim 22 wherein the conditional arrhythmia probabilities 
1 0 are derived using different basic time periods. 

24. The method of claim 1 9 wherein the conditional arrhythmia probabilities 
are based upon past observations of the occurrences of events and arrhythmias 
taken from population data. 

15 

25. The method of claim 19 wherein the conditional arrhythmia probabilities 
are based upon past observations of the occurrences of events and arrhythmias 
taken in real-time from a particular patient. 

20 26. The method of claim 25 wherein the conditional arrhythmia probabilities 
are based initially upon past observations of the occurrences of events and 
arrhythmias taken from population data and each probability is subsequently 
updated from a previous value to a present value with observations taken in real- 
time from a particular patient. 

25 

27. The method of claim 25 further comprising periodically updating a 
conditional arrhythmia probability from a previous value to a present value with 
observations taken in real-time from a particular patient 

30 28. The method of claim 27 wherein a conditional arrhythmia probability is 
updated only if the present value differs by a predetermined amount from the 
previous value. 
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29. The method of claim 28 further comprising testing the amount by which 
the present value differs from the previous value for statistical significance 
before a conditional arrhythmia probability is updated. 

5 30. The method of claim 25 further comprising testing the statistical 

association of a conditional arrhythmia probability with real-time observation 
data. 

3 1 . The method of claim 30 further comprising discontinuing use of a 
1 0 conditional arrhythmia probability in deriving a composite conditional 

probability if the statistical association is below a specified value. 

32. The method of claim 25 further comprising incrementing or 
decrementing a conditional arrhythmia probability by a specific amount after a 

1 5 prediction time period in accordance with whether the arrhythmia occurred or 
not, respectively. 

33. The method of claim 1 further comprising delivering a preventive 
arrhythmia therapy if the estimated arrhythmia probability exceeds a exceeds a 

20 therapy threshold. 

34. The method of claim 33 further comprising selecting a particular therapy 
to be delivered from a group of one or more available therapy modalities. 

25 35. The method of claim 34 wherein the group of available therapy 

modalities includes delivery of one or more of a plurality of pharmacological 
agents. 

36. The method of claim 34 wherein the group of available therapy 
30 modalities includes modalities selected from a group consisting of cardiac pacing 
in a selected pacing mode, delivery of cardioversion/defibrillation shocks, and 
neural stimulation. 



43 



WO 01/24876 PCT/US00/26679 

37. The method of claim 34 wherein the group of available therapy 
modalities includes issuance of a warning signal. 

38. The method of claim 34 wherein the selection of a particular therapy to 
5 be delivered is based upon ascertaining a physiologic state of the patient and 

deciding whether or not a particular available modality of therapy is appropriate 
for delivery to the patient in that physiologic state. 

39. The method of claim 38 wherein the physiologic state ascertained for 

1 0 selecting a therapy modality includes the prediction time period of the estimated 
arrhythmia probability. 

40. The method of claim 38 wherein the physiologic state ascertained for 
selecting a therapy mode includes the particular conditioning events used to 

1 5 make the arrhythmia prediction. 

41 . The method of claim 38 wherein the physiologic state ascertained for 
selecting a therapy mode includes the presence or not of specific conditioning 
events within a specified prior time period. 

20 

42. The method of claim 38 wherein the physiologic state ascertained for 
selecting a therapy mode includes the magnitude of a measured physiologic 
variable. 

25 43, The method of claim 38 wherein the decision as to whether or not to 
deliver a particular mode of preventive therapy takes the form of a matrix 
mapping of a state vector representing a specific physiologic state to a specific 
therapy or therapies considered most appropriate for preventing an arrhythmia in 
that physiologic state, wherein the matrix mapping is performed using a 

30 translation matrix containing information relating to the appropriateness of a 
particular therapy modality for a given physiologic state. 
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44. The method of claim 19 wherein a separate estimated arrhythmia 
probability is computed for each of a plurality of types of arrhythmias. 

45. A processor-readable storage medium having processor-executable 
5 instructions for performing the method recited in claim 1 . 

46. A processor-readable storage medium having processor-executable 
instructions for performing the method recited in claim 19. 

1 0 47. A processor-readable storage medium having processor-executable 
instructions for performing the method recited in claim 33. 

48. A cardiac rhythm management system, including: 

a sensing module for detecting conditioning events statistically 
1 5 associated with occurrences of arrhythmias; and, 

an arrhythmia prediction module for predicting the occurrence of an 
arrhythmia within a specified prediction time period if an estimated arrhythmia 
probability exceeds a specified threshold value, wherein the estimated 
arrhythmia probability is computed from a conditional arrhythmia probability 
20 associated with the conditioning event that is derived from past observations of 
instances in which the conditioning event occurs alone or together with an 
arrhythmia within a specified time period. 

49. The system of claim 48 wherein the sensing module detects conditioning 
25 events selected from a group consisting of a specific morphology of a waveform 

representing electrical activity of the heart, a specific pattern of activation times 
of different areas of the heart as sensed by a plurality of electrodes, a specific 
sequence pattern of heartbeats with respect to time, a value of a measured 
physiological variable, and a statistic based upon a history of occurrences of 
30 conditioning events. 

50. The system of claim 48 wherein the sensing module detects a plurality of 
conditioning events statistically associated with the occurrence of an arrhythmia 
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and further wherein the arrhythmia prediction module compares a composite 
estimated arrhythmia probability with the threshold value in order to predict the 
occurrence of an arrhythmia, the composite arrhythmia probability being a 
combination of the estimated arrhythmia probabilities associated with each 
5 detected conditioning event. 



5 1 . The system of claim 50 wherein the conditional arrhythmia probabilities 
are based upon past observations of the occurrences of events and arrhythmias 
taken from population data. 

10 

52. The system of claim 50 further comprising an adaptive processing 
module for computing conditional arrhythmia probabilities based upon past 
observations of the occurrences of events and arrhythmias taken in real-time 
from a particular patient. 

15 

53. The system of claim 52 wherein the conditional arrhythmia probabilities 
computed by the adaptive processing module are based initially upon past 
observations of the occurrences of events and arrhythmias taken from population 
data and each probability is subsequently updated from a previous value to a 

20 present value with observations taken in real-time from a particular patient. 

54. The system of claim 53 wherein the adaptive processing module tests the 
amount by which the present value differs from the previous value for statistical 
significance before updating a conditional arrhythmia probability. 

25 

55. The system of claim 52 wherein the adaptive processing module tests the 
statistical association of a conditional arrhythmia probability with real-time 
observation data and discontinues use of a conditional arrhythmia probability in 
deriving a composite conditional probability if the statistical association is below 

30 a specified value. 
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56. The system of claim 48 further comprising a therapy module for 
delivering a preventive arrhythmia therapy if the estimated arrhythmia 
probability exceeds a therapy threshold. 

5 57. The system of claim 56 further comprising a preventive therapy control 
module for selecting a particular therapy to be delivered from a group of one or 
more available therapy modalities. 

58. The system of claim 57 wherein the group of available therapy modalities 
10 includes at least one of: delivery of one or more of a plurality of pharmacological 

agents, cardiac pacing in a selected pacing mode, delivery of 
cardioversion/defibrillation shocks, and neural stimulation. 

59. The system of claim 57 wherein the selection of a particular therapy to be 
1 5 delivered by the therapy control module is based upon ascertaining a physiologic 

state of the patient and deciding whether or not a particular available modality of 
therapy is appropriate for delivery to the patient in that physiologic state. 

60. The system of claim 59 wherein the physiologic state ascertained for 
20 selecting a therapy modality by the therapy control module includes at least one 

of: the prediction time period of the estimated arrhythmia probability, the 
particular conditioning events used to make the arrhythmia prediction, the 
presence or not of specific conditioning events within a specified prior time 
period, and the magnitude of a measured physiologic variable. 

25 

6 1 . The system of claim 60 wherein the therapy control module decides 
whether or not to deliver a particular mode of preventive therapy by performing 
a matrix mapping of a list of detected conditioning events to a specific therapy or 
therapies considered most appropriate for preventing an arrhythmia in that 

30 physiologic state, wherein the matrix mapping is performed using a translation 
matrix containing information relating to the appropriateness of a particular 
therapy modality for a given detected conditioning event, and further wherein the 
mapping results in a score for each therapy modality represented in the matrix 
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that is indicative of the appropriateness of each such therapy modality given the 
list of detected conditioning events, the score being evaluated by the therapy 
control module. 

5 62. The system of claim 60 wherein the therapy control module decides 
whether or not to deliver a particular mode of preventive therapy by performing 
matrix mapping of a state vector representing a specific physiologic state to a 
specific therapy or therapies considered most appropriate for preventing an 
arrhythmia in that physiologic state, wherein the matrix mapping is performed 
10 using a therapy matrix containing information relating to the appropriateness of a 
particular therapy modality for a given physiologic state. 
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